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ABSTRACT 



We investigate the physical properties of the interstellar medium (ISM) in the merging pair 
of galaxies known as The Antennae (NGC 4038/39), using the deep coadded ~411 ks Chandra 
ACIS-S data set. The method of analysis and some of the main results from the spectral analysis, 
such as metal abundances and their variations from ^0.2 to ^20-30 times solar, are described 
in Paper I (Baldi et al. submitted). In the present paper we investigate in detail the physics 
of the hot emitting gas, deriving measures for the hot-gas mass (~T0 7 M Q ), cooling times (10 7 - 
10 8 yr), and pressure (3.5 x 10 _11 -2.8 x 10~ 10 dyne cm' 2 ). At least in one of the two nuclei 
(NGC 4038) the hot-gas pressure is significantly higher than the CO pressure, implying that 
shock waves may be driven into the CO clouds. Comparison of the metal abundances with the 
average stellar yields predicted by theoretical models of SN explosions points to SNe of Type 
II as the main contributors of metals to the hot ISM. There is no evidence of any correlation 
between radio-optical star-formation indicators and the measured metal abundances. Although 
due to uncertainties in the average gas density we cannot exclude that mixing may have played 
an important role, the short time required to produce the observed metal masses (< 2 Myr) 
suggests that the correlations are unlikely to have been destroyed by efficient mixing. More 
likely, a significant fraction of SN II ejecta may be in a cool phase, in grains, or escaping in hot 
winds. In each case, any such fraction of the ejecta would remain undetectable with soft X-ray 
observations. 
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Subject headings: galaxies: peculiar — galaxies: individual(NGC4038/39) — galaxies: interac- 
tions — X-rays: galaxies — X-rays: ISM 

1. Introduction 

Massive stars deeply influence the baryonic component of the Universe. Their ionizing radiation and 
the supernovae (SNe) return kinetic energy and metal-enriched gas to the interstellar medium (ISM) from 
which the stars formed (a process usually called "feedback" ) . Feedback exercises an influence not only on the 
gas-phase conditions in the immediate environment of the clusters hosting those massive stars (e.g. McKcc 
1995; Wiseman & Ho 1998; Pudritz & Fiege 2000), but also on the phase structure and energetics of the 
ISM on galactic scales (e.g. McKee & Ostriker 1977; Norman & Ferrara 1996) and on the thermodynamics 
and enrichment of the inter-galactic medium (IGM) on scales of several Mpc (e.g. Chiang, Ryu, & Vishniac 
1988; Heckman 1999; Aguirre et al. 2001). 

The vast range of spatial scales involved is only one of the various difficulties encountered in attempting 
to study feedback. Even when one restricts the discussion to purely mechanical feedback from SNe and 
stellar winds (commonly termed SN feedback), another difficulty is the broad range of complicated gas-phase 
physics, which includes (magneto) hydrodynamic effects such as shocks and turbulence, thermal conduction, 
and non-ionization equilibrium (NEI) emission processes. Yet another complication is that much of the 
energy and metal-enriched material involved is in the hard-to-observe coronal (T > 10 5 K) and hot (T > 10 6 
K) gas phases, requiring the use of space-based far ultra-violet (FUV) and X-ray telescopes. 

Several attempts at putting quantitative constraints on the metal enrichment caused by SN feedback 
in starburst galaxies have been made in the recent past. From a small sample of edge-on starburst galaxies 
observed by ROSAT and ASCA, Weaver, Heckman, & Dahlem (2000) obtained abundances completely 
consistent with solar values and no evidence of super-solar ratios of the a elements with respect to Fe, as 
required by type-II SN feedback (e.g., Gibson, Lowenstein, & Mushotzky 1997). Moreover, they concluded 
that the technique of measuring abundances through X-ray spectral fitting is highly uncertain because of 
ambiguities in the fits (e.g., a degeneracy between the temperature and metallicity) , and strongly dependent 
on the model choice (see also Dahlem, Weaver, & Heckman 1998; Dahlem et al. 2000; Strickland et al. 2002, 
2004). The advent of Chandra (Wcisskopf, O' Dell & van Speybroeck 1996), with its sub-arcsecond angular 
resolution, allowed a better subtraction of the point sources from the galactic diffuse emission, simplifying 
the spectral fitting. However, the low resolution spectra of the Chandra ACIS CCDs (Weisskopf et al. 1996) 
did not allow a dramatic improvement in the accuracy of abundance measures. Notwithstanding the limits 
of ACIS spectral resolution, Martin, Kobulnicky, & Heckman (2002), observing the dwarf starburst galaxy 
NGC 1569, claimed to obtain ratios of a elements to Fe 2-4 times higher than the solar value. They could 
not put constraints on the individual a element abundances because of degeneracies between metallicity 
and spectral normalization. Super-solar a elements to Fe ratios and super-solar abundances of individual 
a elements (e.g., Mg, Si and S) were also reported in other galaxies observed by Chandra, e.g. the face-on 
spirals NGC 6946 (Schlegel, Holt, & Petre 2003) and M 83 (Soria & Wu 2003). 

While X-ray observations of the halos of nearby edge-on disk galaxies may provide the best single probe 
of the action of mechanical feedback on galactic scales, observations of the disks of face-on nearby spiral 
galaxies may provide useful insights into the physics and enrichment of the hot gas. Such observations allow 
a spatially resolved analysis of the hot ISM in the disk that is much less affected by internal absorption than 
in the case of edge-on galaxies. 
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The Antennae are the nearest pair of colliding, relatively face-on spiral galaxies involved in a major 
merger (D = 19 Mpc for H = 75 km s _1 Mpc -1 ). Hence, this system not only allows a study of the 
properties of the ISM relatively unaffected by internal absorption, but it also provides a unique opportunity 
to get the most detailed insight possible into the consequences of a galaxy merger, such as induced star 
formation and its effect on the ISM. 

The presence of an abundant hot ISM in The Antennae was originally suggested by the first Einstein 
observations of this system (Fabbiano & Trinchicri 1983), and has since been confirmed by observations 
with several major X-ray telescopes (ROSAT: Read, Ponman, & Wolstencroft 1995, Fabbiano, Schweizer, 
& Mackie 1997; and ASCA: Sansom ct al. 1996). The first Chandra ACIS observation of The Antennae 
in 1999 gave us for the first time a detailed look at this hot ISM, revealing a complex, diffuse, and soft 
emission component responsible for about half of the detected X-ray photons from the two merging galaxies 
(Fabbiano, Zezas & Murray 2001; Fabbiano et al. 2003, hereafter F03). The spatial resolution of Chandra is 
at least 10 times superior to that of any previous X-ray observatory, which allows us to resolve the emission 
on physical scales of ~75 pc (for D=19 Mpc) and to detect and subtract individual point-like sources (most 
likely X-ray binaries; see Zezas et al. 2002). At the same time, the spectral resolution of ACIS also allows us 
to study for the first time the X-ray spectral properties of the various emission regions, providing important 
additional constraints on their nature. 

While the first Chandra data set demonstrated the richness of the ISM in The Antennae, the number 
of detected photons was insufficient to study its detailed small-scale morphology and spectral properties. 
A deep monitoring observing campaign of The Antennae with Chandra ACIS has produced a detailed and 
rich data set. A first look at this data-set revealed a complex diffuse emission, with signatures of strong line 
emission coming from some of the hot ISM regions, pointing to high a element abundances (Fabbiano et al. 
2004). This paper focusses on the physical properties of the hot gas and its enrichment by SN explosions. 
In a companion paper we describe the data-analysis procedure, the creation of broad-band mapped-color 
images of the diffuse emission, and the performed spectral analysis (Baldi et al. submitted, hereafter Paper 
I), 

This paper is organized as follows. In §2 we briefly review the main results of Paper I. In §3 we derive 
the hot gas parameters, like pressure and density, from the spectral fits. In §4 we discuss the main results 
derived from the abundance measures, focusing especially on the supernova feedback. In §5 we present some 
further insights into the physical properties of the hot ISM, and in §6 we summarize the main results. 

2. Mapped-color images of the diffuse emission and main results of the spectral analysis 

NGC 4038/39 was observed with Chandra ACIS-S seven times during the period between December 
1999 and November 2002, for a total of <~411 ks. The data products were analyzed with the CXC CIAO 
v3.0.1 software and XSPEC package. Full details about the data processing and reduction are given in Paper 
I. 

Paper I presents a mapped-color image of the diffuse emission in The Antennae, representing contribu- 
tions from three different energy bands: 0.3-0.65 keV, 0.65-1.5 keV and 1.5-6 keV. Before combining the 
images from the three bands, all point sources were subtracted and an adaptive gaussian kernel was applied 
to each of the three images. Figure 1 reproduces the mapped-color image for the reader's convenience. Red 
denotes soft X-ray emission (0.3-0.65 keV), green intermediate-energy emission (0.65-1.5 keV), and blue 
hard X-ray emission (1.5-6 keV). A similar procedure was applied to generate a line-strength map for the 
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emission lines from O+Fe+Ne, Mg, and Si in different regions of the hot ISM. Figure 2 shows this second 
map, where red, green, and blue represent O+Fe+Ne emission (0.6-1.16 keV), Mg emission (1.27-1.38 kcV), 
and Si emission (1.75-1.95 keV), respectively. 

Both the mapped-color image of the diffuse emission and the line-strength map guided us in selecting 21 
spectrally similar regions for a proper spectral analysis. These regions are shown in Figure 3. The method 
used for the spectral extraction and analysis is fully described in Paper I. Table 1 resumes the main results 
of the spectral analysis conducted so far. In this table, the goodness of fit (xled over ^ ne number of degrees 
of freedom), absorbing column density Njj, temperature(s) kT, abundances of Neon, Magnesium, Silicon, 
and Iron, and power-law photon index T (if a power-law component is necessary) are summarized together 
with the "best-fit model" (as defined in §5.2 of Paper I) for each region. 

Almost all regions could be fitted with a single thermal component. Only Region 4b required two 
thermal components with different kT; in this region, we could not constrain the low kT value (although we 
found a minimum in the \ 2 statistics at kT = 0.20 keV), and we could put only a lower limit (kT > 0.54 keV) 
to the higher temperature component. Almost all the single-temperature regions required a temperature of 
kT ~ 0.6 keV to fit the data, except for regions 5, 8b, and 14 (kT ~ 0.3 keV). Considering all the regions 
analyzed, the temperatures (both in the single- and two-temperature cases) range from ~ 0.2 to ~ 0.8 
keV (including the errors). Our analysis yielded metal abundances that have acceptable constraints in the 
majority of regions, and led to the most important, and somewhat unexpected, result of our spectral analysis: 
the abundances of Silicon and Magnesium cover a wide range of values ranging from the very low values 
found in regions 2, 6a, 12a, and 12b (Z ~ 0.2Z Q ) to the very high values observed in Region 5 (Z > 20Z Q ). 
The intrinsic absorption throughout the hot ISM is generally low, with Nh ~ 10 20 cm^ 2 typical and often 
consistent with zero. The exceptions are the southern nucleus (regions 8a and 8b) and Region 7, which 
is significantly obscured also in the optical. In the optical and near infrared, Region 7 corresponds to the 
"Overlap Region," where the most active star formation is now occurring (Mirabel et al. 1998; Wilson et al. 
2000; Zhang, Fall, & Whitmore 2001). These three regions show intrinsic absorption of ~ 1 — 2 x 10 21 cm~ 2 . 

3. Hot-gas parameters 

All the calculations performed in this section use the parameters obtained from the best-fit spectral 
models and listed in Table 1. 

Table 2 gives the detailed emission parameters derived from the spectral analysis of each region. The 
observed flux in the 0.3-10 keV band is listed for both the power-law (if present) and the thermal emission 
in columns (2) and (3), respectively. The emission measure (EM — n 2 V) is listed in column (4). For each 
region we present also the intrinsic (absorption-corrected for Njj) luminosity in the 0.3-10 keV band for 
both the power-law (if present) and the thermal emission (columns (5) and (6)). Column (7) gives the ratio 
between the intrinsic luminosities in the 0.3-10 keV band of the power-law and thermal components. In 
the case of Region 4b, modelled with two temperatures, we were unable to determine a confidence range for 
EM because of the lack of constraints on both temperatures (kT\ unconstrained and fcT 2 > 0.54 keV at la). 
Therefore, we give only the values corresponding to the best-fit kT\ and kT 2 (corresponding to 0.20 keV and 
0.60 keV, respectively) 

^From these parameters and those listed in Table 1, we can estimate various properties of the emitting 
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plasma, such as electron density n e (w riff), thermal energy content E t h, cooling time r c , and pressure p. 
These quantities are given in Table 3. The calculation of these parameters depends on the emitting volume 
and on the filling factor 77 of the hot gas. Since we observe regions projected on the plane of the sky, we 
can only measure the "footprint" of each emitting region. To estimate the emitting volumes we adopted the 
following procedure. For all "single" -temperature regions, we assumed a depth of 200K pc for the emitting 
volume, where K is a variable parameter (the numbers given in Table 3 assume K = 1, corresponding to 
the typical depth of a spiral disk). In the case of Region 4b, where we have two temperatures, we assumed 
that the high- and low-temperature components are in pressure equilibrium, and derived the thicknesses for 
the volumes occupied by these two phases from the best-fit kTs, imposing the condition that they add up to 
a total value of 200-ftT pc. The physical area covered by each region scaled to the distance of The Antennae 
is given in column (2). Column (3) lists the thermal gas electron density n e . This quantity depends both 
on the depth-variable parameter K and on the filling factor 77 of the gas and is proportional to r\~ x ^' 1 K^ 1 l' 1 . 
The thermal energy E t h and the cooling time r c (columns (4) and (5), respectively) were estimated following 
Tucker (1975). Both these quantities are proportional to rf^K 1 / 2 . The thermal pressure of the hot gas, 
expressed in dyne cm~ 2 , is listed in column (6). The pressure (p = 2nkT), being directly proportional 
to the density n e , has the same dependence on the variable-depth parameter K and filling factor r\ as the 
density [j]- 1 ! 2 R- 1 ! 2 ). The filling factor 77 is assumed to be unity in our calculations. A filling factor 77 < 1 
for the hot phase of the ISM is plausible from recent results coming from 3D hydrodynamical simulations 
(77 ranging between 0.17 and 0.44 for SN rate between Galactic and 16 times the Galactic value; see, e.g., 
de Avillez & Breitschwerdt 2004 and references therein). However, our parameters have only a slight de- 
pendence on it (oc -q^ 1 / 2 or oc t] 1 / 2 ), hence our estimates are precise within a factor of two or less. As 
described above for the emission measure EM, in Region 4b we were unable to compute a confidence range 
for n e , E^, T c, an d p because of the indeterminacy of the two temperatures. The values of these param- 
eters for Region 4b correspond to the best-fit values of kT\ and kT2- Metz et al. (2004) found somewhat 
larger thermal energies for regions similar to ours. This probably results from their smaller density estimates. 

The calculated cooling times, in the 10 7 -10 8 yr range, are in agreement with the results of F03. The 
pressure is generally of the order of a few 10~ n dyne cm -2 , but reaches higher values (a few 10~ 10 dyne 
cm~ 2 ) in the two nuclei and also in the regions corresponding to the hot-spot regions Rl and R2 of F03. 
These new values are in excess of those found in F03, possibly reflecting the fact that larger areas were 
averaged over in that work. To test this hypothesis we extracted a spectrum for each of these regions using 
the same region files as in F03. The pressures we obtain for all four regions are, indeed, lower than our new 
estimates in the corresponding regions, although still higher than the pressure estimates in F03. F03 argued 
for a possible pressure equilibrium between the hot gas and the cold molecular clouds. Zhu, Seaquist, & 
Kuno (2003) derived from CO measurements a pressure of 4.2 x 10~ n for the northern nucleus (NGC 4038) 
and of 3.1 x 10~ n dyne cm~ 2 for the southern nucleus (NGC 4039). Our observations suggest that the 
hot-gas pressure is a few times higher than the CO estimate for the northern nucleus and almost an order 
of magnitude higher in the southern nucleus. Such a large pressure difference would imply that shock waves 
are being driven into the CO clouds. The hot ISM may therefore in part be responsible for compressing 
and fragmenting the CO clouds, triggering star formation. However, we must bear in mind the large uncer- 
tainties in the estimated pressures, due to the assumptions we had to make on the emitting volume and the 
small-scale properties of the clouds. 

Table 4 presents our estimates of the total mass of hot gas and of the masses of the individual elements 
Ne, Mg, Si, and Fe for each region. The calculation of these individual-element masses was performed using 
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the Anders & Grevesse (1989) compilation of meteoritic abundances. All values in Table 4 are dependent on 
both K and r\ and are proportional to r q 1 / 2 K 1 / 2 . For both the northern and the southern nucleus our total 
mass estimates are a factor of two smaller than the values previously obtained in F03. This discrepancy may 
be due to differences in the best-fit models and to the sizes of the regions used for the spectrum extraction, 
which were larger in the F03 analysis. The total mass of hot gas derived by adding the contributions from 
all individual regions is ~ 1.3 x 10 7 M Q . For the individual elements we obtain total masses of ~ 3.4 x 10 4 
M for Ne, - 1.3 x 10 4 M for Mg, ~ 2.6 x 10 4 M for Si, and - 1.8 x 10 4 M Q for Fe. 

It is worthwile to note that, if we adopt the smaller distance of ~14 Mpc for The Antennae, as derived 
from the tip of the red giant branch by Saviane, Hibbard, & Rich (2004), some of the above physical 
parameters would change. In particular, thermal energies, luminosities, and masses would be smaller by a 
factor of almost 2, while all other listed quantities would remain unvaried. 

4. Abundances 

The unprecedented high angular resolution and sensitivity of Chandra, in conjunction with the long 
integrated-exposure time of our observations, allows us to analyze in detail the properties of the hot ISM in 
The Antennae, down to regions measuring a few hundred parsecs in size. Based on only the first 75 ksec 
observation of 1999 December, F03 were unable to constrain the variety of metal abundances in these regions. 
Indeed, the measure of hot gas elemental abundances in star forming galaxies has always been a challenging 
task (e.g., Strickland & Stevens 1998, 2000; Dahlcm, Weaver & Heckman 1998; Weaver et al. 2000; Dahlem 
ct al. 2000; Strickland et al. 2002; 2004). However, the new, deeper total exposure has allowed us, within the 
limitations of the low spectral resolution of the ACIS CCDs, to place better constraints on the abundances 
than in the past. The method we have used to determine the abundances and their relative confidence 
ranges, described in Paper I, is conservative enough to take into account both statistical and systematic 
biases that may affect the abundance measures in our low-resolution spectra, especially regarding the biases 
deriving from model choice and from the metallicity-temperature degeneracy. In the past, spatially resolved 
measurements of individual elemental abundances were performed on mcgaparscc-sizc regions in the hot 
intra-cluster medium of galaxy clusters (e.g., Mushotzky et al. 1996; Finoguenov, David, & Ponman 2000) 
and, on a much smaller scale, in elliptical galaxies (e.g. Finoguenov et al. 2002). With the new Chandra 
data for The Antennae, we are now trying to explore the "local" chemical enrichment of the hot ISM by 
supernovae within two interacting galaxies. 

As discussed in Paper I, the elemental abundances of Ne-IX, Mg-XI, Si-XIII, and Fe-L that we mea- 
sure in the hot ISM of The Antennae are generally consistent with the stellar abundances measured from 
optical data (~0.5 solar; Fritze-v.Alvensleben 1998), except for a few regions (4b, 5, 7, 8a and 8b) where 
the measured abundances are clearly in excess of this value and, in the case of Region 5, even significantly 
super-solar. Indeed difficulties in measuring the abundances through X-ray spectral fitting of low resolution 
spectra are well documented in the literature (e.g., Weaver et al. 2000). Our very conservative approach, 
however, did not prevent us from locating peculiar regions with bona-fide super-solar abundances. While all 
the other regions have metallicities consistent with solar (within the errors), Region 5 is showing particularly 
enhanced abundances for three out of the four elements we measured {Z^e > 3.9, Zm 3 > 3.8, Zp e > 1.3, 
at 90% confidence level). Such high X-ray derived metallicities of Fe and especially of the a elements are 
quite rare and have not often been found in previous observations of other galaxies. This could be due both 
to the peculiarity of Region 5 and to the fact that constraining abundances through X-ray fitting of low 



-7- 



resolution spectra is always a challenging task. From a Chandra observation of the face-on spiral NGC 6946 
Schlegel et al. (2003) found an abundance of Si= 3.00+};^ for a single-temperature model fit of the diffuse 
emission in this galaxy. Although the errors they quote are at 90%, they estimated this confidence level 
considering only two interesting parameters, a method a lot less conservative than ours. Hints of super-solar 
values of the a elements abundances were also found in the starburst spiral galaxy M 83 by Soria & Wu 
(2002). However, these authors did not fit the abundances, but simply realized that a single-temperature 
thermal model with the a elements abundances fixed at super-solar values and an a to Fe ratio higher than 
solar (e.g., Zjy e = 2.0, Zug = 2.5, Zgi = 3.0, Zp e = 0.9) gives clearly a better fit (xl e d ~ 0.99) than a 
model with all the elements fixed at solar (xled ~ 1-42). All the other measurements derived from Chandra 
ACIS observations of starburst galaxies or ultra-luminous infra-red galaxies (ULIRGs) and presented recently 
in the literature led to abundances consistent with solar (e.g., the ULIRG sample of Huo et al. 2004; the 
edge-on starbursts sample of Strickland et al. 2004) or clearly subsolar (e.g. NGC 253; Strickland et al. 2002). 

The ratios between elemental abundances in each region can be related to the type of supernovae 
generating the metals. If the elemental enrichment results from Type H supernovae, one would expect [Ne/Fe] 
and [Mg/Fe] values approaching 0.3, and [Si/Fe] values approaching 0.5 on average; the corresponding values 
expected in Type la supernova enrichment are dramatically lower (Nagataki & Sato 1998). We used our 
results to explore the origin of the metals in The Antennae. We calculated the [Ne/Fe], [Mg/Fe], and [Si/Fe] 
ratios of the regions listed in Table 1 and plotted them in diagrams displaying the [Ne/Fe] ratio vs. [Mg/Fe] 
ratio (Figure 4a), [Ne/Fe] vs. [Si/Fe] (Figure 4b), and [Mg/Fe] vs. [Si/Fe] (Figure 4c). In these diagrams, 
we also show the theoretical stellar yields expected from SNe la (red symbols) and SNe II (green symbols), 
taken from the compilations by Finoguenov et al. (2000) and Nagataki & Sato (1998). We also indicate 
the regions occupied by the type-II supernova remnants (SNR) Cas A and N132D, using the abundances 
measured by Favata et al. (1997a, b) from BeppoSAX data. To relate our results to other measurements of 
these ratios in the literature we plot also the values relative to the warm Galactic halo (Savage & Sembach 
1996), to the emission both from the disk and from the northern halo of NGC 253 (Strickland et al. 2002), 
to the X-ray diffuse emission observed in NGC 6946 (Schlegel et al. 2003), and to the averaged value in the 
hot halo of a sample of ten nearby edge-on starbursts (Strickland et al. 2004) . 

Although the uncertainties are often large, almost all regions observed in The Antennae show [Ne/Fe] 
and [Mg/Fe] ratios fully consistent with a SN Il-dominated enrichment scenario. However, the [Si/Fe] ratio 
is clearly lower than expected in a SN II enrichment scenario, and in some cases it could be consistent with 
SN la enrichment, although the uncertainties are quite large also for this ratio. The only two exceptions to 
these trends are regions 5 and 7, both located in the upper right corners of the diagrams. These two regions 
have a [Si/Fe] ratio fully consistent with an enrichment stemming entirely from type-II supernovae. Region 
7 is the "Overlap Region" of the two merging galaxies, where the most active star formation is presently 
occurring (Mirabel et al. 1998; Wilson et al. 2000; Zhang et al. 2001). It seems plausible that recent episodes 
of star formation led to the injection of hot enriched gas into the ISM almost exclusively through SNe II. 
Region 5 is the region with the highest abundances of both Mg and Si. Similarly to Region 7, Region 5 
has a high relative star-formation strength and star-formation rate per unit area (see Table 6 and discussion 
below) . 

There are five regions with very low [Si/Fe] values located in the lower part of Figures 4b and c ([Si/Fe] 
< 0, even encompassing the confidence range). These low [Si/Fe] values appear to be inconsistent with SN II 
enrichment and more typical of SN la enrichment. We caution that this surprising result may be affected by 
poor statistics. Typically, in these regions the Si line is detected at less than 2a significance level, in some 
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cases because of the poor statistics at E > 2 keV that may lead to an underestimation of the Si abundance. 
However note that this deficiency may be an independent confirmation of the failure of type II SN yields 
predictions to fit the spectra (e.g. linking together all the a elements using the type II SNe ratio from the 
compilation of average stellar yields from Nagataki & Sato (1998); see Paper I). 

While there may be some uncertainties both in the SN yield models and (even more) in our abundance 
measures, we may also be seeing the effect of different mean ages of star-formation regions and, therefore, 
of different distributions of SN masses in different X-ray emission regions. This observed Si deficiency may 
also be explained if a significant fraction of the Si, larger than for the other elements, has cooled and is 
locked in dust grains. Silicon is one of the most refractory elements, while Neon is one of the least refractory 
elements, as indicated by depletions derived from interstellar absorption lines (see Savage & Sembach 1996 
for a review of this topic). However, Iron is depleted at least as strongly as Silicon, and there is no supporting 
observational or theoretical evidence to claim that Silicon should be more strongly affected than Fe by cooling 
and locking up in dust. 

We note that the abundance ratios measured in the majority of the NGC 4038/39 regions are consistent 
with ratios observed in SNR N132D. On the other hand, the measured ratios do not agree with the values 
observed in Cas A. The different measured ratios in these two SNRs, both originating from a type-II SN 
event, can give us an idea of the complexity and wide variety of input physics involved in yield predictions 
(e.g., convection, reaction rates, method for initiating shock waves; see Gibson et al. 1997). 

Although our measurements are affected by large uncertainties they represent clearly an improvement 
over the previous efforts aimed at measuring the a to Fe ratio in starburst galaxies. The only exception 
comes from the Strickland et al. (2004) sample of edge-on starburst galaxies. This data point, which has 
smaller error bars than our typical error bars, is consistent with the ratios observed in almost all the regions 
of The Antennae. However, note that their smaller error bars derive from the fact that they come from the 
observed hot halo emission in all ten galaxies of the sample. Moreover, in this sample they cannot put a 
stringent upper limit on the [Si/Fe] ratio, a thing that we were able to do in at least some of our regions. 

4.1. X-ray abundances vs. Optical/Radio SF indicators 

To further explore our results, we searched for correlations with age and star formation indicators for 
the different regions of The Antennae. NGC 4038/39 was observed by J/ST/WFPC2 (F336W, F439W, 
F555W, and F814W filters) in 1996 January (Whitmore et al. 1999). There were also VLA observations (in 
configurations BnA, CnB, and B at 6 and 4 cm) performed in 1997 January and 1998 September (Neff & 
Ulvestad 2000). Based on the 7fST/WFPC2 data and VLA observations, we looked for 5 different types of 
objects within each of our regions, counting: (1) the number #i? of very red, extremely young clusters (Zhang 
et al. 2001, 84 in total), (2) the number #i? of bright young blue clusters (Table 1 in Whitmore et al. 1999, 
50 in total), (3) the number #7 of intermediate age clusters (Table 2 in Whitmore et al. 1999, 30 in total, 
of which 5 unpublished), (4) the number fl of flat-spectrum, thermal radio sources (a > —0.40), and (5) 
the number st of steep-spectrum, non-thermal radio sources (a < —0.40). In addition, brightness estimates 
were performed visually on the HST/WFPC2 images, aimed at estimating the integrated luminosity of each 
region in the U, V, I, and Ha passbands on a scale ranging from to 10. 

^From the published data and our visual estimates it was possible to derive — for each region — four 
quantities that might correlate with ISM metal abundances: a mean cluster age, the relative strength of star 
formation, the star formation per unit area, and the number of SNRs per unit area. The mean cluster age 
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in each region was computed as the number-weighted average of the mean ages of very red clusters (4 Myr) 
and young blue clusters (15 Myr); we ignored the intermediate-age (^500 Myr) clusters because their small 
known numbers (0-2 per region) would have added excessive noise. The relative strength of star formation 
was estimated via the ratio (U + Ha)/{V + I), where "relative" is with respect to the older populations 
contributing light to the V and / passbands. The star formation per unit area was calculated as the sum 
of U and Ha estimates (on a scale of 0-10) divided by the estimated area of each region in kpc 2 . The 
number of SNRs per unit area was defined as the number of VLA steep-spectrum radio sources (from Neff 
& Ulvcstad 2000) divided by the estimated area of the region. 

Table 5 summarizes our visual estimates, while Table 6 is a summary of the derived quantities and 
also includes information about the presence of major star-forming complexes named by Rubin, Ford, & 
D'Odorico (1970) in its last two columns. 

Figure 5 displays the relations between the four star-formation indicators listed above and the mean 
abundance of a-elements in each region. The mean abundance was calculated as the straight average of 
the Ne, Mg, and Si abundances, while the error bars were computed through the propagation of errors (at 
lcr). As Figure 5 shows, there are no strong correlations between the different star-formation indicators 
we calculated and the mean abundances in the same regions. The "star-formation per unit area" and the 
number of "SNRs per unit area" are supposed to be among the best indicators of star formation. Yet, these 
two quantities appear uncorrelated with the metal abundances, while a slight correlation may exist between 
the relative strength of star-formation and the abundance of a-elements. Figure 6 compares the same star- 
formation indicators with Fe abundance. No obvious trend between the several quantities is detected in this 
case either, although a larger spread in Fe abundances seems to characterize regions with higher mean stellar 
ages. 

The hot-gas masses derived from abundances in § 3 may yield some clues toward an interpretation of 
this lack of correlations. Assuming the average SN II stellar yields listed in Nagataki & Sato (1998), we 
estimated the number of type-II SNe necessary to produce the observed metal masses in the hot ISM, and 
the time required to produce that number of SNe. If the resulting time is significantly shorter than 100 Myr, 
then the lack of correlations may either be related to the speed at which mixing occurs in the hot ISM 
or indicate that the metals are somehow removed from the hot phase of the ISM. Conversely, if the SN 
production time is >100 Myr, then mean cluster ages may not be relevant because over periods of 100 Myr 
or longer major fractions of the decaying galaxy orbits are traversed and it is unclear whether — even if we 
were able to derive better local mean ages (i.e., including also intermcdiatc-agc clusters) — any relation could 
be established between regional abundances and the past local star-formation history. 

Table 7 lists the results of our calculations for a range of theoretical models for SN II average yields 
(reference given in column 1). Columns (2) - (5) give the numbers of type-II SNe necessary to produce the 
observed masses of Ne, Mg, Si, and Fe, respectively. From these numbers, we estimated the time necessary 
to produce the required number of SNe II, assuming a constant SN rate throughout The Antennae. Neff 
& Ulvestad (2000) derived a SN rate directly associated with the compact radio sources detected in The 
Antennae of 0.03 yr -1 , assuming a typical SNR emits 6-cm radio emission for 30,000 yr and has the luminosity 
of Cas A. However, applying the Condon & Yin (1990) method (based on the ratio of non-thermal radio 
luminosity to the radio-emitting SN rate in our Galaxy) to the total steep-spectrum radio emission of The 
Antennae, and not only to the detected compact sources, they predict a SN rate of 0.26 yr -1 . The resulting 
times required to produce the quantity of Ne, Mg, Si, and Fe observed in The Antennae through SN II 
enrichment were calculated for both of these estimates of the SN rate and are listed in columns (6), (7), (8), 
and (9), respectively. For both rate estimates the times required are short: in the ^3.4-12 Myr range for a 
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SN rate of 0.03 yr _1 and in the ~0.4-1.4 Myr range for 0.26 yr" 1 . Note that the second estimate is more 
reliable because the SN rate was derived taking into account not only the radio emission from the identifiable 
compact sources, but also the overall steep-spectrum radio emission from the whole galaxy. 

However, if the average density of the ISM is lower than our estimates (e.g. if the disk thickness is 
larger than our assumed 200 pc), so that more mass and more metals are needed to produce the observed 
line intensities, then the above times may increase to ~3-4 Myr (or ^25-35 Myr in the lower SN rate 
estimate). Such time spans may be too short for the metals to become mixed throughout the hot ISM. Note 
also that a filling factor of the hot gas lower than 1 would increase the density estimate and make such time 
spans shorter. 

Indeed, 2 Myr is a relevant time threshold for mixing in a hot gas of the temperature we detect, 
corresponding to average velocities of ~ 500 km s _1 . At 500 km s" 1 , gas travels only ~ 1 kpc in 2 Myr. If 
the time required to produce the metals is actually of this order of <2 Myr, it may be unlikely that efficient 
mixing takes place on a scale larger than a few hundred pc, and it is probably insufficient to destroy the 
correlation. Other factors may have destroyed the correlation as well like, e.g., (i) the presence of a significant 
fraction of the SN ejecta in a cool phase or (ii) in dust grains, or (iii) the escape of the ejecta in a galactic 
wind. Another possibility could be that the other metals are invisible in soft X-rays. Indeed, in a Chevalier 
& Clegg (1985) wind model the metals would be expected to be in much hotter (T ~ 10 7 - 5 K) gas. Therefore 
the hot gas observed could be primarily ambient ISM heated, but not particularly enriched, by Type II SNe. 

5. Physical properties of the hot ISM 

The overall characteristics of the hot, diffuse gas in The Antennae, as noted by Read et al. (1995) and 
F03, are (1) a very high pressure, ranging from 10 to 100 times the pressure of the hot component of the ISM 
in the solar neighborhood (Cox & Reynolds 1987), (2) temperatures in the range 0.3 to 0.7 keV, or about 3 
to 7 times the temperature of the Galactic hot ISM, and (3) an abundance of a-process elements generally 
consistent with the ~ 0.5 solar abundance shown by the stellar populations (Fritze-v.Alvensleben 1998). 

There are, of course, uncertainties both in the model fits and in the path length that must be assumed 
to convert emission measures to pressures. In § 3 we assumed a 200 pc thickness for the emitting region and 
a filling factor of 1, unless explicit factors are given. Where two temperatures are required, we assumed the 
components to be in pressure equilibrium and their combined path length to be 200 pc. 

In the spectral analysis above we have pinned down the metallicity parameters for specific regions, and 
we find abundances of Ne, Mg, and Si generally ranging from ~0.2 Zq to Zq, along with a few cases of 
overabundance, represented mainly by five regions: 4b, 5, 7, 8a and 8b. 

Three of these regions can be defined as "Hot Diffuse Emission Regions" (4b, 7, and 8a), requiring fitting 
either with two temperatures, one of which consistently to be above 0.7 keV, or with a single temperature 
above 0.5 keV. On the other hand, the "Cooler Diffuse Emission Regions" (5 and 8b) are both well fit with 
a single temperature around 0.3 keV. 
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5.1. Thermal Energy Supply from type II SNe 

The five regions overabundant in metals listed above represent an heterogeneous subsample of our X-ray 
diffuse emission regions. 

Region 4b lies in the X-ray ridge that runs from NW to SE. Together with regions 4a, 6a and 6b, it is 
bounded to the NW by a strong H I arm and to the SE by a region of strong CO emission. Star formation 
in this region is modest as judged by the number of young clusters (Whitmore et al. 1999), and radio 
continuum emission is correspondingly weak (Hibbard et al. 2001). The mean age of the clusters is about 10 
Myr (Whitmore & Schweizer 1995). Thus this region is relatively quiescent and presents a relatively simple 
case for analysis. 

On the other hand, Region 7 corresponds to the most active burst of star formation currently occuring 
in the Antennae, as observed both in the optical and in the near infrared (Mirabel et al. 1998; Wilson et al. 
2000; Zhang et al. 2001). It also corresponds to the "Overlap Region" between the two merging disks, and 
is discussed in more detail by Fabbiano et al. (2004). 

Regions 8a and 8b are strongly influenced by the presence of the southern galactic nucleus and present, 
together with Region 7, the highest intrinsic absorptions measured in the various spectral fits. 

Finally, Region 5 stands out for its very high abundances and high concentration of young star clusters. 
It, too, is discussed by Fabbiano et al. (2004). Its high abundance of a-elements compared to iron is consistent 
with enrichment by Type II SNe. Its low temperature (see Table 1) could be the result of a relatively high 
local gas density or of relatively effective entrainment of cool ISM gas into the low density SN-heated gas. 
It may evolve into a state similar to that of the neighboring regions 4 and 6 over the course of a few times 
10 7 years. 

Following a similar procedure to the one described in § 4, we derived first the mass of each chemical 
element from the spectral fits and then the number of type II SNe required to produce the observed quantity 
of Ne, Mg, Si, and Fe in the five above regions. The required number of SNe can yield an estimate of the total 
supply of thermal energy available to each region. Our estimates were performed by averaging the numbers 
computed from each of the SN-yield models considered in Table 7 and propagating the errors. Assuming an 
average energy supply of 10 51 erg per SN, we then obtained the thermal energies listed in Table 8. For each 
region the thermal energy supply derived from the observed quantity of Ne, Mg, Si, and Fe are shown in 
columns (2), (3), (4), and (5), respectively. 

If we compare these numbers with the thermal energies derived from the spectral fits (listed in Table 3), 
we can see that the observed thermal energy is systematically lower than the energy supplied by SNe. Note 
that filling factors less than unity would not affect this difference since both thermal energies and masses 
scale as rj 1 ^ 2 . This discrepancy is generally about a factor of 5, with the exceptions of Regions 5 and 8b, 
where the observed thermal energy is lower by about two orders of magnitude. Although the large errors in 
the metallicitics of Region 8b can make this discrepancy less dramatic (and comparable with the other three 
regions), Region 5 docs not present such large uncertainties and clearly stands out with respect to the other 
emission regions. It is likely that in Region 5 most of the thermal energy of the SNe that enriched the hot 
gas (<~ 99%) has been lost either through thermal conduction and radiative cooling or through a wind. 

Radiative cooling of the hot gas can be clearly ruled out by the cooling times listed in Table 3, which are 
definitely too long. On the other hand, the conductive cooling time scale for a volume of size L, temperature 
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T, and density n is 



bnkTL 8 „-5/2 r 2 

<cond = IQ-5T 5 / 2 (T/L)L 2 = X 05 pc 



where To. 5 is the temperature in units of 0.5 keV and L pc is the length scale in parsecs. For n ~ 0.05 cm~ 3 
(the typical electron density we observe in The Antennae X-ray diffuse emission), a temperature of 0.5 keV, 
and a length scale L ~ 30 pc, this formula yields a conductive time scale of ~ 1000 yr, assuming full Spitzer 
conductivity (Spitzer 1962). Thus, conduction may have been able to dissipate a significant fraction of the 
thermal energy produced by Type II SNe in Region 5 and in the other regions we examined, if hot and cool 
gas are mixed on a scale comparable to that in the Galaxy. 

Loss of energy by a wind seems to be a possibility, too. However, a wind is supposed to carry away 
metals as well. Therefore, one would require the energy loss to the wind to be ~ 100 times more efficient 
than the mass loss (at least for Region 5). This would suggest a wind temperature ~ 100 times higher than 
the temperature of the gas we observe, which looks implausible. 



6. Summary and Conclusions 

In this paper and in Paper I we have performed an extensive study of the X-ray properties of the diffuse 
emission of The Antennae (NGC 4038/39), analyzing the entire 411 ks exposure obtained with Chandra 
ACIS-S. Confirming the results of F03, based on the first of the seven observations used in the present 
study, we report a spatially and spectrally complex hot ISM. Thanks to our deep data, we also detect clear, 
spatially variable emission lines of Ne, Mg, and Si, in addition to the Fe-L blend, which indicate super-solar 
metal abundances in a few regions of the hot ISM. In this paper we have tried to investigate the origin of 
these metals. 

To summarize: 

1) We derived physical parameters of the hot gas, such as density, thermal energy, cooling times and 
pressure. This hot gas has cooling times lying in the 10 7 -10 8 yr range, in agreement with the results of F03. 
At least in one of the two nuclei the hot-gas pressure is significantly higher than the CO pressure, implying 
that shock waves may be driven into the CO clouds. Summing the hot-gas masses determined for each region 
yields a total measured mass of ~1.2 x 10 7 M Q for the hot ISM. 

3) Fitting the Fe-L, Ne-IX, Mg-XI and Si-XIII emission, we find significant metal enrichment of the hot 
ISM. Metal abundances are generally consistent with solar, but reach extremes of ~ 20 — 30 solar in Region 
5 and are significantly subsolar in a few regions. The total hot gas masses derived for Fe, Ne, Mg, and Si 
are 1.8 x 10 4 M , 3.4 x 10 4 M , 1.3 x 10 4 M , and 2.6 x 10 4 M , respectively. 

4) Comparison of elemental ratios with those expected from Type la and Type II SNe suggests that 
SNe II are mainly responsible for the metal enrichment of the hot ISM in The Antennae. However, we 
measure [Si/Fe] ratios mostly lower than those predicted for SN II yields. A similar inconsistency was also 
found in Paper I, linking together all the a elements using the type II SNe ratio from the compilation of 
average stellar yields. Although there surely are considerable uncertainties in the SN yields predicted by 
theoretical models, it is also true that our spectra are affected by low statistics at the energy of the Si line 
(E~2 keV) . This suggests that uncertainties in the measurements of the Si abundance may be the primary 
cause of this discrepancy. If our measurements are correct, our result may point to a depletion of Si in the 
hot ISM, perhaps from more efficient cooling of this element. 
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5) We report a remarkable lack of correlations between our abundance measurements and stellar age 
indicators, estimated from cither the Hubble observations of stellar clusters (Whitmore & Schweizer 1995; 
Whitmore et al. 1999) or VLA SNR estimates (Neff & Ulvcstad 2000). The time required to produce the 
observed quantities of metals through type-II SN explosions is probably <2 Myr, implying that efficient 
mixing is unlikely to be the main agent destroying the expected correlation between abundances and radio- 
optical star-formation indicators. This may point toward the presence of a significant fraction of SN II ejecta 
in a cool phase, in dust grains, or escaping in a wind of hot gas, undetectable at soft X-rays wavelengths. 

6) We have analyzed the energy input coming from SN explosions in the 5 regions with metallicitics 
significantly higher than the average stellar abundance (~0.5 solar; Fritze-v.Alvensleben 1998). In all of 
these regions we observe a thermal energy lower than expected from their metallicities and SN enrichment. 
We conclude that the missing energy has been dissipated through processes different than radiative cooling 
(radiative r c ~ 10 7 — 10 8 yr), most likely through conductive cooling. 

We thank the Chandra X-ray Center DS and SDS teams for their efforts in reducing the data and for 
developing the software used in the data reduction (SDP) and analysis (CIAO). We thank D.-W. Kim for 
useful discussions. We also thank the anonymous referee for his careful analysis of the paper and for giving us 
helpful suggestions for the presentation of our results. This work was supported in part by NASA contract 
NAS8-39073 and NASA grants G01-2115X and G02-3135X. F.S. acknowledges partial support from the 
NSF through grant AST 02-05994. 
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Fig. 1. — Mapped-color broad-band smoothed image of The Antennae. The red color channel indicates 
emission from the 0.3-0.65 keV band, green from the 0.65-1.5 keV band, and blue from the 1.5-6 keV band. 



Fig. 2. — Mapped-color line-strength image of The Antennae. Red represents emission from the Fe K-linc 
blend (plus O and Ne), green from the Mg line, and blue from the Si line (for details, see text). The 17 
regions used for a preliminary subdivision of the hot ISM (performed in Paper I) are marked in white. 



Fig. 3. — 0.3-6 keV image of The Antennae, obtained after removing the >3u point sources. The 21 regions 
used for the spectral analysis of the hot ISM are marked in white. 
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Fig. 4. — Abundance-ratio diagrams for regions of the Antennae hot ISM: (a) [Ne/Fe] vs. [Mg/Fe]; (b) 
[Ne/Fe] vs. [Si/Fe]; (c) [Mg/Fe] vs. [Si/Fe]. Abundance units refer to meteoritic abundances of Anders & 
Grevesse (1989). The errors are at 68%, and brackets represent the logarithmic values. Each one of the 
regions in Table 1 is plotted with a black dot, empty if the detection significance of the Mg or Si line in 
the X-ray spectrum is less than 2a. Upper limits are plotted with black arrows. Red symbols correspond 
to ratios expected if all the elements were generated exclusively by SNe la, while green symbols refer to the 
average stellar yield of SNe II. "F00" indicates values adopted in Finoguenov et al. (2000), while the other- 
labels refer to the theoretical work listed in Nagataki & Sato (1998). The magenta and cyan symbols indicate 
abundance ratios for Cas A and N132D, derived from the metallicities observed by BeppoSAX (Favata et 
al. 1997a, 1997b). Blue symbols and lines refer to determinations of abundance ratios in other Galaxies: the 
warm Galactic halo (star symbol or dashed line), NGC 253 (filled square for the northern halo, cross for the 
disk emission), NGC 6946 (lower limit on [Si/Fe] represented as a dash-dotted line) and the Strickland et al. 
(2004) sample of edge-on starbursts (filled circle). 
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Fig. 5. — Diagrams of the mean abundances of a-elements in each region of the Antennae hot ISM plotted 
versus (a) mean age of the star clusters, (b) relative strength of star formation, (c) star formation per unit 
area, and (d) SNRs per unit area. 
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Fig. 6. — Diagrams of the Fe abundance in each region of the Antennae hot ISM plotted versus (a) mean age 
of the star clusters, (b) relative strength of star formation, (c) star formation per unit area, and (d) SNRs 
per unit area. 



Table 1. Best-Fit Parameters (lcr errors for all interesting parameters, see text for details) for the Diffuse Emission of The Antennae. 
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^ r-, ^+ 1 84 




A2 


5 


31.4/35 


9l + il) Bi) 

u - JJ -_o.91 


30+ lu)7 


8.86™ 


19 07+^ ^ 

±y.u( _ 14 q s 


inn oo + ^10.^b 


„ „. +87.43 
^■ Bi -0.95 


1.88*+' 


Bl< + > 


6a 


74.6/61 


74+ 4 tis 


0.66_ 07 


o.i6i u :5 b a 


U " iS -0.28 


o.26ll!; 4 26 


0.27 ± 0.09 


-o.6i^ ; 3 b 9 


Bl 


6b 


95.5/80 


1 47+1 ■'■f 

-1.47 


0.61 ± 0.04 


72 + ' 1 - 44 

u - ' z -0.39 


94+°- 4!J 

u,y -0.37 


o.solU.fg 


43+" 12 

u -*°-0.09 




Al 


7 


86.9/65 


13.77™ 


56+" 12 

U - 00 -0.16 


, or+3.85 


1 93+ 3 0u 

1,ao -l .19 


2.74+ 3 / . 2 


46+"- 43 

u -^ u _0.18 


1.88< + ) 


Bl<+> 


8a 


102.5/90 


8.66I 4 81 


61 + " " 4 


1 50+ lfiS 
1 - OU -0.87 


S9 + (, -88 
U - Oz -0.53 


61 + " Vf< 
u.oi_ 44 


70+ u :l1 

u - '^_o,19 


1.88<+> 


Bl> + ' 


8b 


82.5/71 


24 61+ 11 ™ 

Z *- D1 -11.43 


r\ qo + O.28 

u.oo_ 05 


1 5 57 + °° 
J-0.OI_14.62 


°- y0 -5.33 


12 44 +oc 
"■ M -12.04 


S.lli^f 


3 61 +II B4 

J - Di -0.70 


Bl 


9 


120.6/108 


o gn"l">i-Sti 

Z ' SZ -2.82 


62 + 01 ' 4 




76+ 6B 

u - '°-0.43 


0.841^ 


0.59±«;S 


1.88<+) 


Bl( + » 


10 


128.2/93 


< 2.04 


63+" 04 

U - D,:> -0.05 


33+ 1141 


41+°- 3i 

u '^ 1 -0.28 


u - lu -0.41 


25 + IU)7 

U '^O-0.05 




Al 


11 


57.2/54 




0.62 ± 0.06 


59+ u ' J!> 


, r,y + U.94 

-0.75 


U - 8 '-0.63 


44+"- 17 

u -*^-0.14 




Al 


12a 


55.0/46 


< 2.30 


°- 60 -0.09 


ll + "- 4fi 


35+ :!y 

u - oo -0.32 


20 + " B;> 


„ 97 +o.uy 

U - Z '-0.07 




Al 


12b 


24.5/28 


< 5.29 


0.61 ± 0.06 


< 0.38 


i4+ u<i5 


< 0.82 


n OK + l).12 

°- 26 -o.n 




Al 


13 


32.2/22 


„ „ +ia.u4 

u -° -0.32 


U.OO_ 35 


2.57±^ 2 


, 0O + 1E.-18 

i.oo_ 1 26 


1 54 +13«.7I 


47+^-88 
u -^' -0.30 


1.88 (+) 


Bl*+) 


14 


28.8/21 


O- u <-3.07 


n oy+0.39 
u -° ' -0.08 


unconstr. 


imconsir. 




unconstr. 


1.88<+' 


Bl<+) 


15 


44.9/38 


3 95 + 10 - 92 


59+ 09 


1 09+ 552 

- L - ul '-1.09 


0.89+ 2 ;^ 


U - J '-0.37 


U - ' Z -0.34 


1.88<+) 


Bl(+) 


16 


67.0/78 


n _ c + 4.i)7 

°' 76 -0.76 


62 + " 04 


32 + "- 4B 


52+ :!S 

U ' OZ -0.32 


O.56±" ; 3 4 


34+ U Of! 

u -°^-0.07 




Al 


17 


90.5/73 




62 + " 10 
u.oz_ 34 


96+ 14 74 


„ Q ,+28.48 

y -2.24 


i.ieit;^ 


0.84+ J - 2 2 


1.88*+' 


Bl(+> 



'**XSPEC Models: Al — wabs(wabs(vapec))\ A2 — wabs(wabs(vapec+vapec))\ Bl — wabs(wabs(vapec+powerlaw)); B2 — wabs(wabs(vapec+vapec+powerlaw)). 
^'Photon index frozen. 



Tabic 2. Emission Parameters for the Diffuse Emission of The Antennae. 



.1. 



Region # 
(1) 



F 

(xlO 



10 keV 

" 14 cgs) 



(2) 



rptherra 
^0.3-10 keV 
(XlO" 14 cgs) 

(3) 



em(*~> 

(xlO 61 cm- 3 ) 
(4) 



jp.l. 

-^0.3-10 keV 

(xlO 38 cgs) 
(5) 



r therm 
-^0.3-10 keV 

(XlO 38 Cgs) 
(6) 



rP-i. 

3-10 keV 
i therm 
^0.3-10 fceV 

(7) 



1 






Z 






3 






4a 






4b (+) 






O 


0.44 ±0.21 


6a 


1.68± 32 3 f 


6b 






7 

8a 
8b 
9 


1 
1 

1 


O1+0.42 
zi -0.50 
co+0.46 
D,:> -0.49 

oy + 0.64 

°'-0.40 
q4 +0.61 
94 -0.64 


10 






11 






12a 






12b 






13 
14 
15 







ic+0.35 
1D -0.16 

zo -0.17 

da -0.33 


16 






17 


1 


Ui -0.50 



+0 

+0 

+0 

1.14t° 

+0.32 
0.22 

+0 
-0 
+0. 


+1 

-0 
+2 



o.8o: 

0.65^ 
0.871 
1.141 

0.27±no?, 2 

0.76 
1.92 
3 

0.91_ 
3.151S 

+7. 



2.38 
4.46 
2.86 
1.46 
1.49 
0.74 
0.55 
0.49 
0.81 
2.81 
0.96 



1.65 
1.55 



^+2.49 
-0.74 
hi. 11 
J -1.03 
-+2.02 
-1.31 

4.08+ 3 ' 20 



3.79Z 
2.851 
3.651 



0.97. 
0.85 
8.41 



1.05 

4.17 

1.03 
0.48 
+3.87 
1.43 

+7.20 
-1.61 
+5.7 
-2.5 
f 6.96 
-1.19 
+6.3 
-3.1 



2.76 

10.6 
2.08 

12.4 

11.2tf 7 

4-76l?;i 

1.36 
1.06 
+ 1.36 
-0.70 
+2.35 
1.43 
+0.67 
0.05 
+2.92 
-1.79 
+3.8 
1.7 
+4.32 
1.22 



6.24 
3.11 
1.91 
0.09 
2.17 
10.3 
2.03 



1.93 
7.26 



+0.94 
0.93 
+ 140.12 
5.66 



6.37 
8.12 
6.10 
8.87 



+2.21 
2.63 
+2.29 
-2.44 
+4.51 
2.81 
+2.79 
2.92 



0.68 
1.09 

2.75 
4.64 



1.51 
0.68 
+0.83 
0.74 
1.36 
1.55 



7 +2.3() 



3-47%. 67 

1.10 
0.97 
+2.05 
1.37 
+4.09 
1 33 
1.51 in 0+7.4 
1.01' lu - z -7.0 
+4.13 
1.94 
+4.02 
-1.33 
+5.4 
3.2 
+ 15.98 

3.54 
+ 10.2 

4.3 
+49.4 
-8.5 
+ 10.8 



2.81 
3.78 
5.23 
l.27±\%\, 10 2 
3.41 
8.57 
16.0 
6.10 
18.7 
14 
21.6 



-O.O 

12.41: 2 ; 3 



+2.16 
2.30 



6.81 
6.44 
3.19 
2.40 
2.12 
4.08 
12.5 
4.65 



+5.61 
1.91 
1.42 
1.17 
+ 1.37 
0.72 
+2.87 
-1.80 
15.90 
1.26 
+5.45 
-3.34 
+4.5 
-2.0 
f9.93 
-2.80 



0.566 
0.847 

1.044 
0.434 
0.412 
0.411 



0.283 
0.514 
0.674 

0.998 



MEmission measure defined as EM = n 2 V. 

(+)ln the case of two-temperature gas, the quantities are listed separately for both components: lower and higher temperature, respectively. 



For the calculation of EM we assumed the best-fit value of kTi and kT2\ it was not possible to determine a confidence range for EM because 
of the indetermination in both temperatures 



Table 3. Hot-Gas Parameters. 



Area n e E t ^ t c p 
(kpc 2 ) (xl(T 2 cm- 3 ) (xlO 53 erg) (xlO 7 yr) (xlO^ 11 dyne cm- 2 ) 
(2) (3) (4) (5) (6) 



1 


2.260 


_ _|_ 1 r A 


n 

20-8±|;° 


l9.0t} 4 ;I 


- .-. . -4-4 fi 

10.4l 4 ;8 


2 


0.597 


,-, -1-1 f>1 


^ . -, -1-3 4fi 

9.4lt 3 ;« 




17.9t 7 | 


3 


3.361 




22.0tJ° ° 


18.4± 23 ; 7 


7.43± 3 ; 38 


4a 


0.960 


8.50I 2 ; 8 7 


14.5±£§ 


8.79ll;^ 


. _ . 10 

17 1 ' 


4b(*) 


1.453 


21.2, 7.08 


0.44, 17.0 


1.10, 5.28 


13.6 


5 


0.483 


r ^o + 2 67 


Z - Z4 -1.16 


2.08l?;^ 


r 0^+4 39 


6a 


1.405 


- „ - J-2 1 


26.4+J; 4 


9 - 76 l!'.35 


21-3l3.| 


6b 


2.551 


8 71+i' 3 ,!? 

°- ' -0.95 


38 2+?'? 


7 57+ 4 '° J 9 


17 0+ 3 ' 9 

' — 2.9 


7 


2.351 


4.47l*;g 


16.6± 21 9 7 


8.62+. 38 5 f° 


8.01™ 


8a 


0.856 


14 5+ 3 ' 5 


21 4+ 6 - 8 


3.63t 2 ; 7 4 


28.3+^ 


8b 


0.756 


6.84± 7 ' 43 




1 03+ 8,31 
1 - U >J_0.92 


7.22™ 


9 


4.528 




54.2tif; 5 2 


7 qc+5.97 


13.6+ 4 '! 


10 


7.783 


4 94+0 45 
-0.39 


68.3+^;° 


17 5+ 6 ' 4 
L '-°-5.0 


9.96li;l° 


11 


1.317 


7.84t 2 ;^ 


18.0+. 8 ; 7 


8-3»_4. 86 


15.6+ 7 ' 4 


12a 


0.583 




1 Q Q+2-6 
±J - o -3.0 


6.54t 3 ° 2 


25.9+jj 


12b 


0.266 


14.1+ 2: * 


6 46+ 2 ' 04 


6.42t 4 ;« 


27.5l 8 5 ; 7 7 


13 


0.269 




a oc+3.50 
4.8t)_ 3 . 89 


6.40l 37 8 ™ 


20 4+ 14 ' 8 


14 


0.171 


2 Q8+ 5 ' 70 


1 43+1-74 


2 14+ 9 - 54 
Zil -2.09 


„ co + 17.58 
°-°°-1.77 


15 


0.261 


11 q+6-3 


n ric+5.87 
>J- Z0 _1.59 


9 co+36.50 
z - oz _1.97 


22.5l^ 4 


16 


3.391 


7 20+ 1 ' 21 
'■ zu -0.62 


42 7+ 10 - 4 


IO.8I4; 2 


14 -3 + 23 


17 


5.363 


2 54+ 1 ' 95 


r)o + 25.1 

Z6.»_ 17Q 


1 6 2+ 67 ' 6 


5.04+11? 



(*)ln the case of two-temperature gas, wc calculated n e , E t ^, and t c for both components (assuming the best 
fit temperatures and pressure equilibrium): lower and higher temperature, respectively. It was not possible to 
determine a confidence range for those parameters and for p because of the indetermination in both temperatures. 
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Table 4. Hot-Gas Masses for Hydrogen and Heavy Elements. 



Reg# 


Total Mass 


Ne Mass 


Mg Mass 


Si Mass 


Fe Mass 


(M ) 


(M ) 


(M ) 


(M ) 


(M ) 



1 


5.9 


X 


10 5 


< 1.1 x 10 3 


8.2 


X 


10 1 


< 3.0 x 10 2 


3.1 


X 


10 2 


2 


2.6 


X 


10 5 


1.3 x 10 2 


8.5 


X 


10 1 


5.0 x 10 1 


1.4 


X 


10 2 


3 


7.1 


X 


10 5 


8.5 x 10 2 


3.7 


X 


10 2 


3.6 x 10 2 


4.8 


X 


10 2 


4a 


4.0 


X 


10 5 


3.8 x 10 2 


3.2 


X 


10 2 


9.2 x 10 1 


4.6 


X 


10 2 


4b 


5.3 


X 


10 5 


1.8 x 10 3 


7.1 


X 


10 2 


6.4 x 10 2 


1.2 


X 


10 3 


5 


1.3 


X 


10 5 


2.9 x 10 3 


2.3 


X 


10 3 


1.4 x 10 4 


9.5 


X 


10 2 


6a 


7.0 


X 


10 5 


2.8 x 10 2 


1.8 


X 


10 2 


1.8 x 10 2 


4.9 


X 


10 2 


6b 


1.1 


X 


10 6 


2.0 x 10 3 


9.5 


X 


10 2 


8.7 x 10 2 


1.2 


X 


10 3 


7 


5.2 


X 


10 5 


3.0 x 10 3 


9.2 


X 


10 2 


1.4 x 10 3 


6.2 


X 


10 2 


8a 


6.1 


X 


10 5 


2.3 x 10 3 


4.6 


X 


10 2 


3.7 x 10 2 


1.1 


X 


10 3 


8b 


2.5 


X 


10 5 


9.8 x 10 3 


1.4 


X 


10 3 


3.2 x 10 3 


3.4 


X 


10 3 


9 


1.5 


X 


10 6 


4.8 x 10 3 


1.1 


X 


10 3 


1.3 x 10 3 


2.3 


X 


10 3 


10 


1.9 


X 


10 6 


1.6 x 10 3 


7.2 


X 


10 2 


1.3 x 10 3 


1.2 


X 


10 3 


11 


5.1 


X 


10 5 


7.4 x 10 2 


6.0 


X 


10 2 


4.4 x 10 2 


5.8 


X 


10 2 


12a 


3.9 


X 


10 5 


1.1 x 10 2 


1.3 


X 


10 2 


7.7 x 10 1 


2.7 


X 


10 2 


12b 


1.8 


X 


10 5 


< 1.7 x 10 2 


2.4 


X 


10 1 


< 1.5 x 10 2 


1.3 


X 


10 2 


13 


1.5 


X 


10 5 


9.3 x 10 2 


1.8 


X 


10 2 


2.2 x 10 2 


1.8 


X 


10 2 


14 


2.5 


X 


10 4 


















15 


1.5 


X 


10 5 


4.1 x 10 2 


1.3 


X 


10 2 


5.6 x 10 1 


2.9 


X 


10 2 


16 


1.2 


X 


10 6 


9.5 x 10 2 


5.8 


X 


10 2 


6.7 x 10 2 


1.1 


X 


10 3 


17 


6.7 


X 


10 5 


1.6 x 10 3 


1.8 


X 


10 3 


7.7 x 10 2 


1.5 


X 


10 3 



Table 5. Mean Optical and Radio Observed Properties of the 17 Chandra Regions in NGC 4038/39. 





HST 


HST 


HST 


VT A 

J 1 


VT A 

V 1 jr\ 


Phot 


Phot 


Phot 


Phot 


Rc g # 


WFPC2 


WFPC2 


WFPC2 


HST 


HST 


HST 


HST 




#R 


#B 


#1 


U 


V 


I 


Ha 




1 

1 


n 
U 


U 


u 


1 

I 


1 

1 


Z 


n 
O 


1 
1 


Z 


1 


Q 

o 


U 


r 
.0 





y 


O 


A 


1-7 

I 


3 


2 











2 





2 


A 


1 


A 


i 


Z 




U 


U 


i 


Q 

o 





Z 


c 



Z 


Q 
O 


U 


o 
z 


l 







Q 
O 


Q 

o 


6 


5 


6 


2 





3 


4 


2 


3 


4 


7 


16 








3 


11 





1 


3 


4 


8 


2 


1 





1 


2 


3 


6 


10 


10 


9 





7 





3 


4 


3 


4 


5 


10 


10 








1 





1 











1 


11 


4 


3 


2 


1 


5 


6 


3 


3 


8 


12 


2 


6 


1 


1 


4 


10 


6 


5 


10 


13 





3 





2 





6 


5 


5 


8 


14 


1 











1 


1 


6 


10 


5 


15 





2 








2 


3 


3 


8 


6 


16 


2 


2 





1 


2 


2 


3 


5 


3 


17 


10 


2 


1 


3 


15 





2 


3 


1 



( a ) Values for Region 1 based only on eastern part (western covered by WFPC2 to < 50%) 



Table 6. Mean Optical and Radio Derived Properties of the 17 Chandra Regions in NGC 4038/39. 





Mean 


Relative 


SF per 


SNRs 


Unit 








Ago 


Strength 


Unit 


per 


Area 


Gal 


Landmark Objects 




(Myr) 


SF 


Area 


Area 


(kpc 2 ) 






l(o) 


4.0 


0.4 


1.10 


0.55 


1.81 


38 


North in "W Loop" 


2 


12.2 


1.8 


27.6 


8.6 


0.58 


38 


"G" complex 


3 


4.0 


0.2 


0.31 


0.63 


3.2 


38 


Nearly pure old disk 


4 


11.3 


0.3 


1.14 


0.0 


2.64 


38 


"H" & outer bulge, dust lane 




10 6 


2 3 


29 2 


2 1 


48 


38 


Z / O l j tllSK 


6 


10.0 


1.6 


2.1 


0.8 


3.8 


38 


2/3 "E," strong SF in NE disk 


7 


4.0 


1.0 


1.7 


4.6 


2.4 




''Overlap region" 


8 


7.7 


0.8 


8.9 


1.4 


1.46 


39 


"A," nucleus & bulge of 39 


9 


15.0 


1.4 


3.2 


1.0 


4.0 


39 


"B" + "AA," outer bulge of 39 


10 






0.13 


0.13 


7.8 




Intcrdisk region, ~cmpty 


11 


8.7 


2.3 


10.8 


3.8 


1.3 


38 


"M," south in "W Loop" 


12 


12.2 


1.8 


29. 


5.7 


0.70 


38 


"S" + "T," in "W Loop" 


13 


15.0 


1.4 


56. 


0.0 


0.25 


38 


"R," in "W Loop" 


14 


4.0 


0.4 


32. 


5.3 


0.19 


38 


"J," N part of inner bulge 


15 


15.0 


0.8 


39. 


8.7 


0.23 


38 


"K," S part of inner bulge Sz nucl. of 38 


16 


9.5 


0.6 


1.6 


0.65 


3.1 


38 


"L," outer bulge of 38 


17 


5.8 


0.2 


0.22 


3.3 


4.5 




"Overlap region," most dusty part 



^Values for Region 1 based only on eastern part (western covered by WFPC2 to < 50%) 



Table 7. Number of SN II and Time Required to Produce the Observed Metal Masses. 



Reference ^ 

(i) 


iVsNii(Ne) 


ATsNii(Mg) 


JVsNii(Si) 


JVs N ii(Fe) 


*snii(No) <+) 


tSNIl(Mg) < + ' 


tSNii(Si) (+) 


tSNIl(Fo) ( + ) 


(xlO 5 ) 


(xlO 5 ) 


(xlO 5 ) 


(xlO 5 ) 


(Myr) 


{Myr) 


{Myr) 


{Myr) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


A96 


3.4 


2.4 




2.5 


11.3, 1.3 


8.0, 0.9 




8.4, 1.0 


T95 


1.5 


1.1 


2.0 


1.5 


4.9, 0.6 


3.7, 0.4 


6.6, 0.8 


4.9, 0.6 


T95+M92 








1.5 








4.9, 0.6 


W95;A;1O" 4 Z 


3.6 


3.6 


2.5 


2.5 


12.0, 1.4 


12.0, 1.4 


8.4, 1.0 


8.2, 0.9 


W95;B;1O~ 4 Z 


1.5 


2.0 


2.2 


2.1 


5.1, 0.6 


6.6, 0.8 


7.4, 0.9 


7.0, 0.8 


W95;A;Z 


1.9 


2.0 


2.1 


1.6 


6.3, 0.7 


6.7, 0.8 


7.0, 0.8 


5.3, 0.6 


W95;B;Z 


1.3 


1.4 


1.8 


1.3 


4.3, 0.5 


4.6, 0.5 


6.1, 0.7 


4.2, 0.5 


N98S1 


1.5 


1.1 


1.9 


2.1 


5.0, 0.6 


3.5, 0.4 


6.3, 0.7 


6.9, 0.8 


N98A1 


1.5 


1.0 


2.5 


2.1 


4.9, 0.6 


3.4, 0.4 


8.2, 1.0 


6.9, 0.8 


N98A3 


1.5 


1.1 


2.7 


2.0 


4.9, 0.6 


3.5, 0.4 


8.9, 1.0 


6.8, 0.8 



<*>A96, Arnett 1996; T95, Tsujimoto ct al. 1995; M92, Maeder 1992; W95, Woosley & Weaver 1995; N98, Nagataki ct al. 1998. N98S1, N98A1, and N98A3 
are the models of the spherical explosion and the axisymmetric explosions, respectively. The degree of deviation from spherical symmetry is larger in N98A3 
than in N98A1. 

^'Computed for a SN rate of 0.03 yr _1 and of 0.26 yr — 1 , respectively (see text). 
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Table 8. Thermal energy supplies for the five regions presenting overabundances of metals. 



Rcg# 
(1) 



Bth(Ne) 
(erg) 
(2) 



Bth(Mg) 
(erg) 
(3) 



3* (Si) 
(erg) 

(4) 



^th(Fe) 
(erg) 
(5) 



4b 
5 
7 

8a 
8b 



1.02±^o x 10 55 
1.64±g-|| x 10 55 
1.72t? : 2J x 10 55 
1.30tS:S x 10 55 
5-61±g°27 x !0 55 



9T2™ 4 x 10 54 



3.06 



+51.99 
-2.49 



x 10 



55 



1.23j££ x 10 55 
6.16±I |? x 10 54 
1.87+f 6 ° 4 x 10 55 



5.41^-10 x 1Q54 

1.20±f;|o x 10 56 
1.19±o?| x 10 55 
2.50± 4 ;° 4 x 10 54 
2.66±^ 59 x 10 55 



1.271^8 x 10 55 

1.02±g 2 3 | 4 x 10 55 
6.6ll 6 2 ;g x 10 54 



1.18 



-0.54 
0.33 



X 10 



55 



3.6lt 46 1 ° 1 x 10 55 
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